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ABSTRACT 

We study gravitational instabilities in disks, with special attention to the most massive clumps 
that form because they are expected to be the progenitors of globular- type clusters. The max- 
imum unstable mass is set by rotation and depends only on the surface density and orbital 
frequency of the disk. We propose that the formation of massive clusters is related to this 
largest scale in galaxies not stabilized by rotation. Using data from the literature, we predict 
that globular-like clusters can form in nuclear starburst disks and protogalactic disks but not in 
typical spiral galaxies, in agreement with observations. 

Subject headings: galaxies: disk instabilities - globular clusters: formation - star formation: general 



1. Introduction 

The study of instabilities in disks has a long his- 
tory, following the seminal work of Toomre (1964), 
with a considerable literature on many aspects of 
it. However, relatively little attention has been 
given to the most massive agglomerations that can 
form by the fragmentation of galactic gas disks. 
These most massive agglomerations are of interest 
because they may be the precursors of the most 
massive star clusters known, the globular clusters. 

Globular clusters were until recently viewed 
as exclusively old objects, and cluster formation 
models were therefore based on ideas about early 
stages of galaxy formation. This view changed 
with the realization that elliptical galaxies often 
contain two populations of globular clusters that 
appear to have different origins, suggested to be 
a 'primordial' population and a 'merger' popu- 
lation (Ashman & Zepf 1992). The existence of 
a merger population was confirmed by the Hub- 
ble Space Telescope with the discovery of mas- 



sive young clusters in merging galaxies (Ashman 
& Zepf 1998; Schweizer 1998). A modern theory 
of globular cluster formation should therefore ad- 
dress why globular cluster formation is common 
in some environments such as merging and high- 
redshift galaxies but not in others such as the 
present Milky Way. 

Another question to be addressed is why some 
regions of galaxies are more favorable for the for- 
mation of massive clusters than others. For ex- 
ample, the inner part of the Milky Way galaxy 
contains young clusters with masses up to sev- 
eral times 10 4 M Q , including some quite close to 
the central black hole (Krabbe et al. 1995; Figer 
2008), that are an order of magnitude more mas- 
sive than the typical open clusters found elsewhere 
in the galaxy. 

The Giant Molecular Clouds (GMCs) in which 
open clusters currently form in the outer Milky 
Way have masses of order 1O 6 M0, but they pro- 
duce clusters with masses of only several times 
10 3 M Q . The formation of globular clusters with 
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masses of 10 5 — 10 6 M@ therefore requires either 
much more massive super-GMCs (Harris & Pu- 
dritz 1994) or a much higher star formation effi- 
ciency. In reality, a combination of these effects is 
probably involved. 

Several scenarios have been proposed for the 
formation of massive super-GMCs, involving for 
example a hot primordial plasma (Fall & Rees 
1985), collisions between normal GMCs (Harris & 
Pudritz 1994), or confinement by high pressures 
in mergers (Ashman & Zepf 2001). Relatively lit- 
tle attention has been given to the possible role of 
disks in the formation of globular clusters, with 
the exception of Larson (1988, 1996), but this 
possibility now seems worth further investigation 
given the evidence for massive cluster formation 
in rotationally supported disks in merging galax- 
ies (Ashman & Zepf 1998). 

In this Letter we study the gravitational insta- 
bility of galactic gas disks, with special attention 
to the most massive clumps that can form in them. 
We also consider which galaxies or environments 
will produce the largest unstable clumps and are 
therefore most favorable for massive cluster forma- 
tion. We start by reviewing the theory of instabil- 
ities in disks in §2, and continue with its applica- 
tion to the formation of globular-type clusters in 
§3. We discuss the most promising environments 
for globular-cluster formation in §3.1, and in §4 we 
summarize the main implications of our work. 

2. GRAVITATIONAL INSTABILITY IN 
DISKS 

For illustration we consider the simplest possi- 
ble case: a uniformly rotating thin sheet or disk. 
Linear stability analysis of such a system yields 
the dispersion relation for small perturbations in 
surface density (Binney & Tremaine 1987) 

uj 2 = 4n 2 - 2ttG£|/c| + k 2 C 2 , (1) 

where C s = \J~jj^ is the sound speed, S is the sur- 
face density, and Q is the angular rotation speed. 
To see the implications of equation [I] we first con- 
sider the limit of a non-rotating sheet with f2 = 0. 
In this case the sheet is unstable (oj 2 < 0) for wave- 
lengths A = 27r/fc larger than the Jeans length, 
A jeans = C 2 /GYi. Thus pressure can only stabi- 
lize the sheet at small scales. Another limiting 



case is a pressureless sheet with C s = 0. Equa- 
tion [T] then shows that only perturbations with 
A < A rot = 7r 2 G£/f2 2 are unstable, so that rota- 
tion can only stabilize the sheet at large scales. 
The existence of a maximum length scale not sta- 
bilized by rotation was first derived and discussed 
by Toomre (1964). 

Clearly neither pressure nor rotation can by it- 
self stabilize the sheet, and there is a range of un- 
stable length scales limited on small scales by ther- 
mal pressure (at the Jeans length Aj ean s) and on 
large scales by rotation (at the critical length set 
by rotation, A rot ). All intermediate length scales 
are unstable, and the most rapidly growing mode 
has a wavelength 2Aj cans . Only a combination of 
pressure and rotation can stabilize the sheet, and 
this happens when the range of unstable wave- 
lengths shrinks to zero; according to equation [1] 
this occurs when Aj ea ns > A ro t/4, which is equiv- 
alent to C s r2/G£o > 7r /2, the classical Toomre 
(1964) criterion for the stability of a uniformly ro- 
tating disk. 

Such stability analyses apply only to fluid disks 
that are well described by a simple equation of 
state (EOS), such as an isothermal or a poly- 
tropic EOS. Unfortunately, the real interstellar 
medium in galaxies is highly complex and not 
well described by a simple equation of state. The 
complex dynamics and thermodynamics of the 
real interstellar medium produce large fluctuations 
and sharp transitions in its properties, as hap- 
pens for example in molecular clouds, that change 
the Jeans mass by orders of magnitude and are 
not well described by a simple EOS. A medium 
with complex structure and dynamics may there- 
fore have no well-defined Jeans length, and there 
may be no real lower limit on the sizes of the 
self-gravitating structures that can form until the 
much smaller thermal Jeans scale is reached in 
molecular cloud cores. This thermal Jeans scale 
may be the next smaller scale that has any clear 
physical basis and where a relatively simple EOS 
may again apply; stability analyses give length and 
mass scales in molecular cloud cores of around 0.1 
pc and 1 Mq. 

The maximum length scale set by rotation A ro t 
is however still a meaningful quantity even when 
the ISM is not well described by a simple EOS, 
because A ro t depends only on the surface density 
and angular frequency of the disk and not on its 
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small-scale physics. If a disk has a well-defined av- 
erage surface density, it should have a well-defined 
maximum length scale for instabilities. This rota- 
tional maximum scale may then be the most rel- 
evant length scale for the large-scale dynamics of 
a disk with complex small-scale physics. In the 
Milky Way, A rot is typically of the order of 1 kpc. 
At intermediate scales between 0.1 pc and 1 kpc 
there may not be any preferred scale, and it has 
been suggested by Larson (1979, 1981) and oth- 
ers (Ballesteros-Paredes et al. 2007 and references 
therein) that the structure and dynamics of the 
ISM on these intermediate scales may be roughly 
self-similar and described by power laws, as in a 
turbulent cascade. 

Since the maximum unstable length scale set by 
rotation is the only clearly meaningful length scale 
for a gas disk with complex physics, we focus on 
the characteristic mass associated with this scale. 
The resulting maximum unstable mass, calculated 



as M™ ax 
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For typical parameters in the Milky Way disk, 
equation [2] gives a maximum unstable mass of the 
order of a few times 1O 6 M , similar to the scale 
of giant cloud complexes in our Galaxy. Equation 
[2] gives a mass scale similar to that derived previ- 
ously by many authors (Balbus 1988; Elmegreen 
1994, 2002; Kim & Ostriker 2001) on the as- 
sumption that it is the Jeans mass -Mj cans = 
£ g as(Ajeans/2) 2 that determines the masses of gi- 
ant cloud complexes, where the Jeans mass is 
calculated assuming that on large enough scales 
the ISM can be modeled as an isothermal fluid 
in which random cloud motions take the place of 
thermal motions. This mass scale is again about 
10 6 — 10 7 M Q , and it is similar to the rotational 
maximum mass because the analysis assumes that 
disks are marginally stable (Q ~ 1), implying that 
the Jeans scale and the rotational scale are similar. 

A significant radial variation is implied by [2j 
since both E gas and fl are functions of radius 
in galactic disks. For disks with flat rotation 
curves that have density profiles steeper than 
£ gas oc R~ 4 / 3 , the maximum clump mass de- 
creases with radius. For example in the Milky 
Way, the maximum mass decreases with radius 
from about 10 6 Mr at the solar radius to 3 x 



10 4 M Q at 20 kpc from the galactic center. 

In the more general case of non-uniform rota- 
tion, equations [T] and [2] are easily generalized by 
replacing 2f2 by the epicyclic frequency k, which is 
given by k 2 (R) = R^m- + 4S1 2 . Since the epicyclic 
frequency n varies only between f2 and 2Q for cen- 
trally concentrated disks, we will continue to use £1 
because it has a more simple and intuitive mean- 
ing. For a review of results of more general gravi- 
tational stability analyses see Larson (1985). 

The effect of rotation on the properties of the 
massive clumps that form in disks is illustrated 
by the simulations of nuclear gas disks by Escala 
(2007). These simulations study the evolution of a 
gas disk with a mass of 3 x 10 8 Mq and a diameter 
of 600 pc, varying only the constant rotation veloc- 
ity of the disk and studying its effect on the disk's 
evolution. This was done by varying a fixed exter- 
nal potential in such a way that the disk is always 
rotationally supported, but with a range of rota- 
tion velocities. Here we compare the development 
of the instabilities in two disks with rotational ve- 
locities v rot of 110 and 321 kms -1 . Further details 
of the simulations are given by Escala (2007). 

Figure Q] illustrates a representative stage in the 
evolution of the two disks and shows the density 
distribution in the disk midplane after 0.4 orbital 
periods. The disk with a rotational velocity v mt 
of 110 kms -1 is shown in Figure [1^ and the disk 
with u rot = 321 kms -1 in Figure [TJd. The plot- 
ted region in the x-y plane is 0.5 kpc on a side. 
Both disks develop a complex clumpy structure, 
and the figures illustrate how rotation affects the 
growing instabilities. In the run with lower ro- 
tation (Fig. [1^) the clumps are denser and more 
numerous, while in the more rapidly rotating disk 
(Fig. [TJd) the clumps are less marked and filaments 
are more prominent, reflecting the disruption of 
some clumps by rotation. This is an expected re- 
sult since Figure [TJd has a smaller M^ ax , so that 
rotation can here disrupt clumps that are able to 
survive in Figure [T^l; this explains why the clumps 
are denser and more numerous in the run with less 
rotation. 

The gas clumps in Figure QJl have typical 
masses of several times 10 6 M©, in rough agree- 
ment with the maximum clump mass of 1.5 x 
10 7 Mq predicted by equation [2] For the disk 
with higher rotation (Fig. Q]b), the clumps have 
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smaller typical masses of several times 10 5 M@ , 
also in agreement with the value of 4 x 10 5 Mq 
predicted by equation [2] These values are valid 
for the early evolution of the disk before much of 
its gas has been consumed by star formation. At 
later times, as E gas decreases, less massive clumps 
are formed. 

3. Applications: The largest proto-globular 
clusters 

The condition for rotational support of a homo- 
geneous uniformly rotating gas sheet against the 
gravity produced by the total mass enclosed within 
any radius R, M tot = M Btar + M D m + M gas , can 
be written 

n2 = ^ p _ (3) 

where rj = M gas /M tot is the ratio of the gas mass 
to the total enclosed mass. For a uniformly ro- 
tating, rotationally supported gas sheet, the max- 
imum clump mass predicted by equation [2] can be 
written 

Mci i - 3 10 M& Wmq~ W • 

(4) 

Equation 0] shows that the formation of massive 
clusters requires not only a large gas content, as 
in the Milky Way, but also a large gas fraction 
(large rf) ; the formation of massive clusters is thus 
strongly favored in regions where the gas consti- 
tutes a large fraction of the mass. 

In the local universe this condition is satisfied 
in ULIRGs, which are typically characterized by 
an ongoing burst of star formation. Recent ob- 
servations have shown that massive clusters are 
in fact currently forming in these systems (Ash- 
man & Zepf 1998). In these starburst systems, 
not all of the star formation is confined to a sin- 
gle well-defined disk, especially during the early 
stages of a galaxy merger before it settles down 
into an ordered flow, as in the Antennae system. 
However, the maximum mass derived above is due 
to a purely local effect and not dependent on large- 
scale ordered motion. The mass scale is set by 
rotation via the local balance between centrifugal 
and self-gravitational forces, and it is not strongly 
sensitive to the geometry; equation [2] should there- 
fore remain approximately valid within geometri- 
cal factors. For a typical nuclear starburst disk 



with a mass of a few times 10 9 Mq and a gas frac- 
tion -q of 1/5 (Downes & Solomon 1998), we pre- 
dict the formation of clumps with masses up to 
10 8 Mq. The formation efficiency of massive clus- 
ters is very uncertain, but even with an efficiency 
of only 1% we predict the formation of clusters 
with masses up to 1O 6 M0, characteristic of the 
most massive globular clusters. 

In normal disk galaxies like the Milky Way, 
the total gas mass is similar but the gas frac- 
tion is only about 1/30 of the total mass. We 
then predict clumps with masses of order a few 
times 10 6 M Q , as mentioned in §2, that are not 
massive enough to form globular clusters. How- 
ever, conditions for massive cluster formation be- 
come much more favorable in the protogalactic gas 
disks that are expected to exist during the early 
stages of disk galaxy formation, where the gas is 
expected to constitute 1/5 or more of the total 
mass (77 ~ 0.2). Such protogalactic disks, like for 
example the clump cluster and chain galaxies ob- 
served at high redshift (Elmegreen & Elmegreen 
2005), would then also be very favorable sites for 
massive cluster formation. 

Thus we see from equation 0] that there are two 
favorable environments for the formation of mas- 
sive clusters: (1) nuclear gas disks in starburst 
systems, and (2) protogalactic disks. This expec- 
tation agrees with the scenario for globular clus- 
ter formation proposed by Ashman & Zepf (1992), 
whereby globular clusters form either in the early 
stages of galaxy formation or in galaxy mergers 
and interactions. This scenario has been success- 
fully tested against several observed properties of 
globular cluster systems, such as the typically bi- 
modal age and metallicity distributions of globu- 
lar clusters in elliptical galaxies and the presence 
of young globular clusters in merger remnants. 
Moreover it has been recently found that the old- 
est dwarf galaxies have the highest globular cluster 
frequencies (Peng et al. 2008), which is expected 
since systems that formed earlier would probably 
have had higher gas fractions and therefore they 
were more favorable sites for cluster formation. 

3.1. Environments for Massive Cluster 
Formation 

The formulation presented in this Letter also al- 
lows us to predict which present environments are 
most favorable for forming massive clusters, given 
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the total gas mass and the gas fraction in galax- 
ies, quantities that can easily be extracted from 
observations of many systems. Figure [2] shows 
the total gas mass of a sample of galaxies plotted 
against the gas fraction r\ = M gas /M to t, where the 
black dots are nuclear starburst disks, the open 
circles are spirals, and the triangles are dwarf ir- 
regular galaxies. The data for spiral galaxies are 
from Kent (1987), Giraud (1998), and Heifer et 
al. (2003). For dwarf irregular galaxies the data 
are from Cote et al. (2000), Van den Bosch et al. 
(2001), and for starburst galaxies they are from 
Downes and Solomon (1998). The parallel dotted 
lines show the maximum clump mass predicted by 
equation 31 with values varying from 10 4 Mq to 
10 8 A/q; the solid line is for a maximum clump 
mass of 10 7 M Q . 

We expect globular clusters to form in clumps 
with masses above 10 7 M , so that galaxies above 
the solid line in Figure [5] should be good candi- 
dates for massive cluster formation. If massive 
clusters form with an efficiency of 1% in GMCs 
(Lada & Lada 2003), this corresponds to cluster 
masses of 10 5 M or more, comparable to globu- 
lar cluster masses. Figure [5] shows that with one 
exception, only galaxies classified as nuclear star- 
bursts lie on or above the 10 7 M Q line, in agree- 
ment with the fact that only starburst galaxies 
show ongoing massive cluster formation. The one 
exception is the spiral galaxy NGC 2403, which 
has many giant HII regions and is forming mas- 
sive stellar clusters like starbursts (Drissen et al 
1999), in agreement with our prediction. 

Figure [2] also shows that dwarf irregular galax- 
ies should not form massive clusters even though 
most of their baryonic mass is in gas. This is be- 
cause the more important quantity is the gas frac- 
tion, which is only between 1% and 10% in these 
systems because they are dark-matter dominated. 
In spiral galaxies, although the total amount of gas 
is comparable to that in starburst galaxies, the gas 
is more widely distributed and the gas fraction in 
the star-forming region is therefore smaller, so we 
do not expect massive cluster formation in most 
spirals, in agreement with observations. 

The position of the Milky Way in Figure is 
denoted by MW and the position of the Large 
Magellanic Cloud is denoted by LMC. Both the 
Milky Way and the LMC are predicted to have 
maximum clump masses of several times 10 6 Mq, 



which for a cluster formation efficiency of 1% is 
consistent with the fact that they are both cur- 
rently forming clusters with masses up to several 
times 10 4 M Q (see Figer 2008 for the Milky Way 
and Larson 1988 for LMC). However, given the un- 
certainty in the cluster formation efficiency, what 
is more significant is the prediction that the Milky 
Way and the LMC should form clusters two orders 
of magnitude less massive than the most massive 
clusters that form in starbursts, as is observed. 

4. DISCUSSION 

In this paper we have studied gravitational in- 
stabilities in disks, with special attention to the 
most massive clumps that form because they are 
expected to be the progenitors of globular-type 
clusters. The maximum unstable mass is set by ro- 
tation and depends only on the surface density and 
orbital frequency of the disk, unlike other mass 
scales such as the Jeans mass that depend on the 
complex gas physics. This mass scale is due to 
a purely local instability and not dependent on 
large-scale ordered motion. The maximum unsta- 
ble mass is therefore a well-defined quantity even 
if the interstellar medium is not well described by 
a simple equation of state. 

The maximum clump mass can be expressed in 
terms of the total gas mass and the gas fraction in 
a galaxy, and this formulation makes it clear that 
environments with a high gas fraction are the most 
promising places to form massive clusters. Using 
data from the literature, we predict that massive 
globular-like clusters can form in nuclear starburst 
disks and protogalactic disks but not in typical 
spiral galaxies, in agreement with observations. 

The scenario proposed here relates massive 
clusters to the largest scale in galaxies not sta- 
bilized by rotation, which is the only scale in- 
termediate between stars and galaxies that has 
a clear physical basis. There is no well-defined 
'Jeans mass' on these intermediate scales, and 
the next smaller scale that has any clear physi- 
cal basis is the thermal Jeans scale in molecular 
clouds, which is related to the masses of individual 
stars. The next larger physical scale is that of the 
galaxy itself, so we predict three physically well- 
motivated scales corresponding to stars, massive 
stellar clusters, and galaxies. 

Another application of studies of the stabil- 
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ity of nuclear gas disks, not discussed here, is 
to their likely role in feeding supermassive cen- 
tral black holes in galaxies (Escala 2007). This 
requires the outward transfer of angular momen- 
tum, and if gravity is the most important force 
involved, the same mass concentrations that form 
massive clusters may also play an important role 
in the outward transfer of angular momentum and 
the growth of a central black hole. We plan to in- 
vestigate this problem further in ongoing work. 
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Fig. 1. — a) Density distribution in the plane of the gas disk, in units of 1.19 xlO~ 21 g cm -3 and coded on a 
logarithmic scale, at time t = 0.4 t or b for the disk with rotational velocity of 110 kms" 1 . The plotted region 
in the x-y plane is 0.5 kpc x 0.5 kpc. The figure shows a complicated multiphase structure, characterized 
by high density clumps that are embedded in a less dense medium, b) Same as in a) for the disk with 
rotational velocity of 321 kms -1 . The figure also shows a complicated multiphase structure as in a), but 
now characterized by filaments that are embedded in a less dense medium 
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Fig. 2. — Total gas mass in a galaxy plotted against the gas fraction M ga ,s/Mtot- The black dots are 
nuclear starburst disks, open circles are spirals and triangles are dwarf irregulars. The parallel dotted lines 
indicate predicted clump masses from equation Q] of lO 4 Af , 10 5 M©, 10 6 M Q and 10 s Mq. The solid line 
is for 10 7 M Q . MW denotes the position of the Milky Way in the plot, and LMC the position of the Large 
Magellanic Cloud. 
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